Fermionic ground state at unitarity and Haldane Exclusion Statistics 
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We consider a few-particle system of trapped neutral fermionic atoms at ultra-low temperatures, 
with the attractive interaction tuned to Feshbach resonance. We calculate the energies and the 
spatial densities of the few-body systems using a generalisation of the extended Thomas-Fermi 
(ETF) method, and assuming the particles obey the Haldane- Wu fractional exclusion statistics 
(FES) at unitarity. This method is different from the scaled ETF version given by Chang and 
Bertsch (Phys. Rev. A76,021603(R) (2007)). Our semiclassical FES results are consistent with the 
Monte-Carlo calculations of the above authors, but can hardly be distinguished from their over all 
scaling of the ETF result at unitarity. 
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I. INTRODUCTION 

There has been a lot of interest in a dilute gas of neutral 
fermionic atoms at ultra-cold temperatures both exper- 
imentally |l| and theoretically [2|. In general, the low- 
energy properties of the gas are determined by the scat- 
tering length a, the number density n, and the tempera- 
ture T of the gas ( the effective range ro is small, so that 
ro/|a| — > as a becomes large ). When the attractive 
interaction between the atoms is increased continuously 
by magnetic tuning from weak to strong, the scattering 
length a goes from a small negative to a small positive 
value. In between, there is a zero-energy two-body bound 
state, and \a\ is infinite. The gas is said to be at unitar- 
ity in this situation, and the length scale a drops out. 
The behaviour of the gas is expected to be universal at 
unitarity Q • Understanding the ground state of the sys- 
tem in this limit is a challenge to many-body theorists 
as originally discussed by Bertsch 0]. Experimentally, if 
the temperature is small enough, a BCS superfluid is ob- 
served at the weak end, and a BEC condensate of dimers 
at the strong end This was predicted long back by 
Leggett Q, who extended the BCS formalism in a novel 
fashion to analyse the physical situation. The BCS to 
BEC transition is found to be smooth, with no disconti- 
nuity in properties across the unitary point. 

Chang and Bertsch Q have recently presented an ab 
initio Green's Function Monte-Carlo (GFMC) calcula- 
tion of the energy and density of N = 2 — 22 trapped 
fermionic atoms in a harmonic potential. The atoms are 
interacting via a short range central two-body potential, 



with its strength adjusted to yield a zero-energy two- 
body bound state in free space. The many-body proper- 
ties of this system are expected to be independent of the 
shape of the two-body interaction in such a set-up. 

In an earlier paper Q we obtained the energy per par- 
ticle and the chemical potential of a noninteracting gas 
of atoms at finite temperatures obeying fractional exclu- 
sion statistics (FES) We assumed that at unitarity, 
the effect of the interaction could be simulated by FES 
for the bulk properties of the system. Our results, with 
the choice of one free parameter in FES, were found to 
be in good agreement with theoretical MC calculations 
for a free gas 0, [13] and the experimental results in a 
trap 
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12j . Since the number of atoms, N, was taken 
to be large, no finite-N corrections were needed in our 
semiclassical calculations. This is not the case in the 
present paper, where N is taken to be small. The pur- 
pose of this paper is to test whether the FES hypothesis 
gives improved results when the finite-N corrections are 
incorporated in our calculations. 

We briefly recall the rationale for using FES. As is 
well known, FES is realised by the Calogero-Sutherland 
model [13j in one dimension [14| . In two dimensions, the 
kinetic and potential energy densities of fcrmions inter- 
acting with a zero-range potential scale as the square of 
the spatial density, and obey FES in the mean-field ap- 
proximation 15|. A hint that Haldane statistics is also 
realised for cold atoms in three dimensions at unitar- 
ity comes from the observation that the total energy per 
atom of the gas may be obtained by scaling the kinetic 
energy term by a constant factor [2|. This is however 
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a necessary, but not a sufficient condition. A further 
hint comes from the fact that the second virial coefficient 
of the gas at unitarity is temperature independent [r^ . 
In exclusion statistics, the scale-invariant interaction be- 
tween atoms alters the ideal Fermi (Bose) values of the 
(exchange) second virial coefficient +(— )2 -5 / 2 by adding 
an interacting part |16j | . When FES is incorporated in the 
T = Thomas-Fermi (TF) method, it gives the same 
expression as the scaled density functional approach of 
Papenbrock 2l| . This constitutes the bulk of the smooth 
part of the energy. It is the small next order term in the 
Extended Thomas-Fermi (ETF) expression that behaves 
differently with scaling or FES. The former was done by 
Chang and Bertsch [5]. In this paper, we incorporate 
FES in an improved version of ETF to test if it can dif- 
ferentiate between the two alternatives when compared 
with the GFMC results. 

In the next section, we first summarise the TF and 
ETF results. Both these have limitations at the classical 
turning point, where the spatial density behaves discon- 
tinuously. To rectify this, we make use of a modified semi- 
classical method [18( that gives a continuous variation of 
the density across the turning point. Our semiclassical 
results incorporating FES are next compared with the 
GFMC calculations of Chang and Bertsch Q for fermions 
trapped in a three dimensional oscillator potential. We 
find that FES results are consistent with the many-body 
GFMC results. 



II. SEMICLASSICAL CALCULATIONS 
INCORPORATING FES 

There has been much interest amongst theorists to 
calculate the properties of a gas in the unitary regime 
(fc/|a| >> 1), where kj — (S^n) 1 / 3 is the Fermi wave 
number of the noninteracting gas. This is a challenging 
task, since there is no small expansion parameter, and a 
perturbative calculation cannot be done. In particular, at 
T = 0, the energy per particle of the gas is calculated to 
be H = £ | -i^jj- , where £ is a constant scaling factor [2j . 
A MC calculation gives £ ~ 0.44 [19]. The experimental 
value is about 0.5, but with large error bars [201 ]. A sim- 
ilar relation may be obtained assuming a noninteracting 
gas obeying FES, with a statistical occupancy factor g. 
At T = 0, FES gives 
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where we have included a spin degeneracy factor of 2 
and the factor 1 j g is the FES occupancy factor at T = 
with g = 1 as the fermionic limit. The modified Fermi 
momentum kf, from above, is kt = g x ^ 3 kf, where kf is 
the fermi momentum of the noninteracting Fermi gas. It 
also follows that the energy per particle of the unitary gas 
is given by jj — g 2 ! 3 | -^f . Comparing with the scaled 
version, we see that the scaling factor £ in a Fermi gas 
is related to the statistical parameter g by the relation 
£ = g 2 ^ 3 '. In a three-dimensional isotropic harmonic trap, 



a similar scaling of the TF expression gives 2l| , in units 
of Hu>, 
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In FES, an identical relation is obtained, with £}/ 2 re- 
placed by g 1 / 3 - The scaled TF spatial density is also 
identical to the FES expression when this replacement is 
made : 

I r 2 \ 3/2 

2~P 



Ptf{t) 



1 



(2) 



where I = y/ (h/muS). The above expression is valid for 
r < r*o, where tq = V^lfigN) 1 / 6 is the classical turning 
point. For r > ro, the TF density is zero. To implement 
finite-N corrections, one has to consider ETF [22j . Chang 
and Bertsch [5| scale the energy expression for ETF by 
the same over all factor as in TF (denoted by ETF'), 
where as FES yields a different expression [3] (in units 
of Hui): 
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Although ETF gives a reasonable description of the 
smooth part of the energy, it fails to do so for the spa- 
tial density. In fact, the ETF density diverges at the 
turning point. To give a consistent description of both 
the energy and the spatial density, we adopt a method 
where a selective summation of the higher ord er g radient 
terms of the Wigner-Kirkwood series is made [18| . For a 
harmonic trapping potential V(r), retaining terms up to 
third order in /3, the Bloch density C (r, f3) incorporating 
FES is given by 

1 ,2m. 
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The spatial density is obtained by taking the inverse 
Laplace transform of C(r, (3)/ ' (3 with respect to the chem- 
ical potential /i, which we denote by p = C~ 1 [C(r 1 /?)//?]. 
It is the cubic term in the exponent that makes the den- 
sity continuous across the classical turning point. Simi- 
larly, the energy is given by E = p,N — £~ 1 C(r, /?)/ /3 2 . 
The inverse Laplace transformations are carried out by 
the saddle-point method. The quality of the approxima- 
tion is tested by applying the method to N noninteract- 
ing spin-1/2 fermions (g = 1) in a harmonic potential. 
The result for the energy is plotted as a function of N 
is plotted in Fig.l. To facilitate the comparison, the TF 
energy is subtracted out from the quantum as well as 
the semiclassical results for Eetf and E. Note that the 
shell effects in the energy as well as the density are not 
reproduced in the semiclassical calculations. 
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FIG. 1: Plot of the energy, E — Etf vs particle number for 
fermions (i.e., g — 1). The red (solid) curve corresponds to 
the exact calculation in a harmonic oscillator while the green 
(dashed) and blue (dotted) curves correspond the calculations 
based on ETF and resummation methods. 

In Fig. 2, we compare the GFMC results [H] of the en- 
ergy for N — 2 — 22 atoms with the various semiclassical 
calculations. For the latter, the scaling factor in Eq.Q 
is taken to be £ = 0.48, that corresponds to g — 1/3 for 
the ETF Eq.gl), and also for E. Our choice of £ = 0.48 
is very close to that of [5| , and corresponds to a a not 
too different from the value of 0.29 chosen earlier |6j . It 
is seen from Fig. 2 that all the semiclassical methods fare 
well, and it is not possible to distinguish the scaled results 
from the FES ones. A clearer comparison is made in the 
inset of Fig. 2, where the large TF term given by Eq.JT]) 
is subtracted out from the energies. Even then, it is not 
possible to assert the relative superiority of over all seal- 
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FIG. 2: Plot of the energy vs particle number with g = 1/3. 
The data points refer to the GFMC calculation of Chang and 
Bertschji]. The pink(dotted) curve corresponds to the en- 
ergy calculation using the resummation method. The green 
(dashed) and blue(short-dashed) curves correspond to the 
standard ETF calculation with FES (see Eq.@) and the 
ETF' approximation given by Chang and Bertsch (see eq.©). 
Inset shows the energy after subtracting the TF contribution 
as in Fig.l but choosing g — 1/3. 

ing to FES. We suggest that a distinction may possibly be 
made if a larger range of N values are spanned by a MC 
many-body calculation. An interesting aspect of GFMC 
results (see inset in Fig. 2) is the odd-even oscillations in 
energy. In Fig. 3 the calculated density for N = 20 parti- 
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FIG. 3: Plot of radial density for TV = 20 as a function of 
the scaled distance x. Shown are the density calculated using 
resummation method (red-solid) and the TF density (green- 
dashed) with gr = 1/3. The density of fermions in a harmonic 
oscillator is also shown for comparison (blue-dotted). 
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cles is plotted using the resummation method and is com- 
pared with the TF density. Although there is not much 
to choose between the ETF and the resummation results 
for the energy, the density in the latter case is distinctly 
superior and appears to agree with the smoothed part of 
the density calculated by the the GFMC method. It also 
reproduces the tail beyond the turning point, which is 
not possible in TF or ETF approximations. 

To summarise, We have considered a few particle sys- 
tem of trapped and interacting neutral fermionic atoms 
at ultra-low temperatures. The energy and spatial den- 
sity of this system is calculated semiclassically assuming 
the particles obey the Haldane-Wu fractional exclusion 



statistics (FES) at unitarity. The semiclassical FES re- 
sults are consistent with the Monte-Carlo calculations of 
Chang and Bertsch [5(, but can hardly be distinguished 
from the over all scaling of the noninteracting energy that 
is commonly used at unitarity. However, it is interesting 
to note that both at finite temperature Q and at zero 
temperature the FES frame work yields reasonably good 
results. 
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